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Abstract 

Expressions for the absorption spectrum of a nucleus in a three- and a two-dimensional crystal 
respectively are obtained analytically at zero and at finite temperature respectively. It is found 
that for finite temperature in two dimensions the Mossbauer effect vanishes but is replaced by 
what we call a Quasi-Mossbauer effect. Possibilities to identify two-dimensional elastic behavior 
are discussed. 
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I. INTRODUCTION 



In nuclear reactions, the atoms usually are treated as free, i.e., unbound. The chemical 
bonds between each other are neglected. This is justified with the argument that the typical 
energy in a nuclear process is generally much higher than the chemical binding energies (1 
to 10 eV) and moreso than that of lattice vibrations (10~^ to 10~^ eV)^. In return, most of 
the detailed qualities of the nucleus are of secondary importance in solid state physics. 

These assumptions are not correct without limitations. Consider for example the absorp- 
tion or emission of 7-photons by nuclei. During absorption processes (emission analogously) 
the nucleus first is in a ground state and makes a transistion into an excited state. The 
energy difference between the ground and the excited state is not constant. It has a certain 
distribution which is linked to the life-time of the excited state by Heisenberg's uncertainty 
relation. For these reactions the difference has values between 10~^^ to 10~^ eV— . Thus, it 
is considerably smaller than the energy of chemical bonds. 

During the absortion event, the atom would recoil due to the conservation of the sums 
of the momenta of the atom before the reaction and the absorbed photon. The thermal 
movement of an atom in a gaseous absorber leads to a Doppler broadening of the absorption 
line and the recoil to a shift of its maximum. Typical recoil-energies lie in the range from 
10~^ to 10"*"^ eV— . This interval also comprises the energy of chemical bonds, phonons, and 
even energetically lower areas. 

For atoms bound in a system and photon momenta below the threshold of the destruction 
of the chemical bonds, there are different aspects coming into play. For molecular systems 
additional rotational and vibrational degrees of freedom can be important. Simultaneously, 
the Doppler broadening and the recoil of the system decrease due to the enhanced mass of 
the system. 

Making the transition to a crystal, the Doppler broadening and the recoil of the entire 
system become negligible. The entire recoil-energy goes into lattice vibrations. Further, 
there is the possibility that the momentum of the photon is taken up by the entire crystal 
without producing any phonons. Such an event is called recoil-free. In a rigid crystal, the 
absorption spectrum would just show the pure absorption spectrum of the nucleus, because 
no phonons can be produced. For photons with very low energy, i.e., less than 1.5 x 10"*"^ 
eV the recoil-energy is below the energy of the phonons. Classically, no vibrational modes 
of the crystal could be excited, thus the energy and momentum conservation could not be 
satisfied and the related processes could not take place. However, observations show that 
they do take place and that in those situations a part of the events are recoil-free. Energy 
and momentum conservation are guaranteed when averaging over many events^. The effect 
and the fraction of recoil-free events have been named after its discoverer R. L. Mossbauer— . 

The above discussion also plays a role in the scattering of x-rays or neutrons from atoms. 
There, the distinction between elastic and inelastic events is also known. Mossbauer used 
also results by Lamb about neutron scattering^. This is the reason for the name Mossbauer- 
Lamb fraction. 

Of interest is not only the total intensity not contained in the Mossbauer peak but 
also its distribution onto the remaining spectrum. Here, Visscher-^- derived an analytical 
result for the observable spectrum at zero temperature in a three dimensional crystal. He 
also performed numerical calculations for finite temperatures. Systems with almost two- 
dimensional magnetic interactions are known (interaction strength 10^-times stronger in 
the plane than between the planes). Now, it is interesting to identify systems with low- 
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dimensional elastic interactions. Here, one could think about ultrathin films sputtered onto 
the surface of a solid, a liquid, or the boundary layer between two liquids. Thus a mechanical 
decoupling to the substrate should be achieveable at least in directions parallel to the surface. 

In the following section we discuss the theoretical description of these phenomena by 
deriving and solving a differential equation for the spectrum. Subsequently, in section III, 
we evaluate the solution in harmonic approximation for the isotropic Debye model in three 
and two dimensions at zero and finite temperature. In two dimensions at finite temperature 
we find that there are no recoil-free events. The Mossbauer effect ceases to exist in the strict 
sense. However there is still a central peak made up of almost recoil-free events, a Quasi- 
Mossbauer effect. Ways how to identify two dimensional elastic behavior are discussed. In 
the final section we summarise our paper. 

II. THEORY 

Starting from the absorption spectrum of a nucleus inside a rigid crystal lattice ^^^^(ti;^), 
we write the observed spectrum in an elastic lattice S{uj^) as a convolution of the first 
with a shift spectrum: S = I * Srig. Here, is the circular frequency of the converted 
photon. The shift spectrum I{uj) comprises all the information on the modifications of 
the observed spectrum S{uy) due to the dynamics of the lattice. Up to an overall factor 
I{uj) is the spectrum one would observe, if the nucleus' line width was exactly zero. A 
normalisation to one J dujl{uj) = 1 guarantees the conservation of the integrated intensity: 
J duiS{ijj) = J duSrigiuj). 

In a crystal with N Bravais lattice sites, A atoms per elementary cell, and where the 
atoms are elongated in d' dimensions there are M = NAd' vibrational modes with circular 
frequencies VLj with J G {1,2,..., M}. In general, the dimensionality of the crystal d is the 
same as that of the elongations of the atoms d' which does not exclude the possibility of 
doing calculations with d ^ d' . 

The state of the crystal before (j) and after (j) the emission or absorption process can 
be described by a set of M numbers each {ra}*'-^ = {nf : J G {1,2, ...,M}}. The set of all 
possible changes of the state is constructed by combining every possible initial state with 
every possible final state (direct product). The different changes contribute to the shift 
spectrum I{uj) with different weights yl({n}* — > {n}^). The shift spectrum can be expressed 
as the weighted mean over all changes conserving the total energy: 



Due to the linearity of the above expression, the energy difference in equation is 
determined exclusively by the change 



(1) 



{nY {n}! 

In a harmonic crystal the total energy is given by: 




(2) 
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mj = rij — rij 
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of the occupational numbers in the different modes but not from the state {n}* of the 
crystal before the absorption or emission process. The weight yl({n}* {"n}^) is equal to 
the absolute probability for the transition from the state {n}* to the state {n}'^. It is given 
by the product of the probability A({n}*) of the crystal to be in a state {n}* before the 
process and the conditional probability to reach a final state {n}-^ coming from 

an initial state {n}*: 

A{{ny^{ny) = A{{ny)A{„,y{{ny). (4) 
Substitution of equations and (jD) into equation ^ leads to: 



{m} \ J J {ny 

where the elements of {n}-^ are to be replaced in accordance with equation 0. For a 
harmonic crystal in thermal equilibrium before the elementary process, the terms belonging 
to a certain mode are independent from the configuration of the rest of the system. After 
defining the thermally averaged - i.e., averaged with the thermal weight P^i for having rfj 

phonons in mode J before the process^ - probability for the production of mj phonons in 
mode J: 



nij / , jUj inj^irij+mj^ 



where A'^^^^.^ is the probability for a change of the phonon number in mode J from 
ifij to iUj + mji^, the shift spectrum can be reexpressed as: 



{mj} \ J J J' 

The probabilities are normalised to one: Mnj — 1- Consider the shift spectrum 

I^-^^Lu) for only J active modes. The shift spectrum I^'^^cu) for a nucleus in an entirely rigid 
crystal lattice is given by a Dirac-5-distribution: 5{uj). The shift spectrum with all 

modes in action is identical to the actual shift spectrum I{uj). From equation ((7j) we get a 
recursion equation for the partial shift spectra: 



I^'\u:) = Y,^ij^'''\^-^M- (8) 

mj 

An analogous equation is valid for the changes AJ^'^^ = I^'^\uj) — I^'^~^'>{uj) in connection 
with the shift amplitudes a;^^ = A,;^^ — 5mj,o'- 



A/(^)M = (9) 

mj 
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The shift amphtudes a;^^ differ from the probabihties A:^^ only for mj = 0. There the 
respective shift amplitude is smaller by one. For example A:^^ = for all mj 7^ and 
Aq = 1 signifies that no phonons are created or annihilated in mode J; so does a;^^ = 
for all mj. Hence Aq is a meassure for remaining in the previous state and Oq indicates 
the relative change. Because of the additional Kronecker tensor 6ij and the normalisation 
of the probabilities A'^^, the sum over all shift amplitudes amounts to zero: Xlmj '^mj ~ ^■ 
Writing the shift spectrum I^-^^uj) as a sum over all its changes and replacing the changes 
AI^-^^uj) of the shift spectrum by the recursion Q yields: 



J'=l rrij/ 

In case, several active branches are present, the label of the mode J is composed of a 
branch index j and a Bloch vector k: J = {j,k). In a transition from a discrete to a 
continuous description the summation over all modes J is replaced by a summation over 
the different modal branches and an integration over fc-space. The dimensional integral 
is substituted by a one-dimensional integration over the frequencies Q of the modes. As a 
correction, the density of modes i]j{Q) for the corresponding branch has to be taken into 
account. In the following, all investigations will be carried out for a single active branch. 
So, with the volume V of the crystal: 

E - / tSj ^ / (11) 



(27r 

With the dispersion relation Q{k) the density could be given by^: 



J 

6. 



vm = V 1 6[Q - Qik)]0^^. (12) 

It can be interpreted as the total number of modes in the infinitessimally small interval 
+ dQ]. The modes are counted in the entire crystal and not only in a unit volume. 
So, the density of modes grows proportionally to the number of Bravais lattice sites. It 
is to vanish for frequencies outside the intervall [0, f^A/]- When connected to the discrete 
description this must also be the interval in which lie all the modal frequencies Qj with 
J G {1, 2, M}. In the Debye model, Qm is equal to the Debye frequency Qd- For simple 
models, the alternative to marking the shift spectrum I^-^^uj) and the shift amplitudes 
a:^j with the mode index J is to express them as functions of the modal frequency Qj: 
I^-^\uj) = and a;^^ = am{^j)- During the transition from the discrete to the 

continuous description the discrete parameter flj is replaced by the continuous Q. Values 
Sit Q = Qj for all J G {1,2,...,M} are preserved: Qj, J G {1,2,...,M} ^ Q e [0,Q], 
I{Qj,uj) —>■ I{Q,uj), and a^i^j) — > am{^)- Carrying out the continuum limit in equation 
(fTn|l leads to: 

^0 
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Here - for the sake of simplicity - we anticipate a result that will be derived later on: The 
shift amplitudes cim,(^) for fixed argument are proportional to integer negative powers of 
the number of Bravais lattice sites: am(^) ~ iV~l™'l for all m 7^ and cio(^) ~ A^~^. The 
shift amplitudes occur only combined in a product with the density of modes. As already 
stated above, the density of modes is proportional to the number of Bravais-lattice sites. 
So, only the three terms with |m| < 1 survive taking the continuum limit. Two of the 
three shift amplitudes that are left can be expressed as functions of the third: They can 
be eliminated with the help of the normalisation of the shift amplitudes and because we 
are assuming thermal equilibrium: a_m(f^) = e'"^^^^ a+rn{^) , respectively, where (3 = ksT. 
Partial differentiation of equation ()13|) with respect to the maximum frequency VL together 
with the above facts leads to a linear, first-order, non-local, partial differential equation: 



with the starting condition /(0,a;) = 6{uj). A Fourier transformation from the variable lo 
to the variable and a subsequent division through the Fourier transformed shift spectrum 
I{fl, 0) leads to an ordinary first-order differential equation: 



— ln{/(^], 0)} = r/(fi)a+l(^]){e-'3^^(e+^^^ - 1) + (e^^^ - 1)}. (15) 

The actual shift spectrum /(ti;)_is equal to I{Q = Qm, i^)- Taking into account the Fourier 
transformed boundary condition /(0 = 0) = 1, we find by direct integration: 

i{(f)) = exp |y dQ7]{Q)a+i{Q) [e-/^^^(e+*^^ - 1) + (e"*^^ - 1)] | • (16) 

Given the existence of the resulting integrals, the integral in the exponent can be seperated 
into a 0-dependent and a 0- independent part. A Fourier retransformation yields the shift 
spectrum I{uj) in momentum representation: 

J(cu) = e'^ J ^e'<^" exp | dnri{n)a+^ (fi) [e(+^^-/5^)^ + e"'*^] | . (17) 

The factor e"^ with F = J^^' dQr]{Q)a+i{Q){e~'^^^ + 1} is equal to the Mdssbauer- 
Lamb fraction. The exponential series is uniformly convergent on every bounded area of the 
complex plane. So, as long as the magnitude of the exponent is finite, the exponential series 
can be integrated term by term: 



I(^uj) = e~^J2^ I —e'^^i I "f/fiT/(fi)a+i(r])[e(+*^-'3^)^ + e-^<^^]| . (18) 

Reinterpretation of the inner integration as an additional Fourier transformation and 
decomposition of the sum into single terms Iy{uj) results in: 
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1J\ ATI 

where denotes a Fourier transformation and where (27r)^^jF0^_^ is its inverse. 

The evaluation of the first two terms is independent of the form of the density of modes or 
the shift amphtude: 



h{Lo) = e-^ {r]{-Lj)a+i{-ij)e+^^'' + r]{+u;)a+i{+u;)} . (20) 

Terms with u > 2 could be reexpressed with the help of the binomial theorem and by 
noting that products become convolutions after a Fourier transformation. Then one sees 
that the term with index u contains the effect of events on the shift spectrum during which 
u phonons are converted, i.e., created or annihilated. 



III. APPLICATION 

Now, we would like to evaluate explicitely the shift spectrum I{uj) in harmonic approxi- 
mation for the isotropic Debye model. The Debye model is based on the assumption that the 
crystal is a (i-dimensional hypercubic lattice. So, it is a Bravais-crystal in which no optical 
modes exist. Isotropy leads to degenerate modes on all d' = d branches. The frequencies 
follow a linear dispersion relation: VL{k) = sk with the speed of sound s and the absolute 
value k of the wave-number vector k. Additionally, the integration over the first Brillouin 
zone is replaced by an integration over the volume of a sphere in /c-space. The size of the 
sphere is chosen in such a way that it contains allowed wave-number vectors where N 
is the number of ions in the crystal. Evaluation of equation (jl2j) under these assumptions 
leads to: 



r^{n) = d'— i^—j 9inn - n)ein), (21) 

with the Debye frequency fin = sk^. As mentioned above, the Debye frequency replaces 
the maximal frequency Qm- The transition probabilities which are the thermally averaged 
absolute squares of the transition matrix-elements {{n'j + Tn'j}\e^'^^\{n'j}) are to be evaluated 
in a harmonic system. Here, q stands for the momentum transfer and x for the spatial 
coordinate along its direction. After some calculations one finds: 



At 



—p"^ coth 



sinh 



-I -1 



(22) 



with the parameter: 
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and where Im stands for the modified Bessel function of the first kind and Mjo„ for the 
mass of the ion whose nucleus absorbs the 7-photon. For small arguments, we could use 
the lowest order approximation to the Bessel functions: lm{w) = (iy/2)™/m! + 
However, the argument of the Bessel function is divergent for a vanishing circular frequency 
Q. Due to the periodic boundary conditions there exists a minimal circular frequency Qi 
which, for a constant volume density of atoms n, depends on the total number of atoms 
according to: ^li = SN^^^'^. The constant of proportionality S depends on the geometry 
of the lattice and for a (i-dimensional hypercubic lattice is given by S = 27csn~^^^. For 
large N, the argument of the Bessel function becomes: w = 2g^iV*^^~'^)/'^/(Mjo„/3S^). In 
three dimensions it goes to zero as goes to infinity. In two dimensions it is independent 
of the number of atoms A^, so for the validity of the lowest order approximation we have 
to postulate additionally: 2g^/(Mjo„/5S^) << 1. In one dimension the argument w is not 
boundend from above and so the lowest order approximation is not adequate. In every case 
where the lowest order approximation is justified, we have: 



^+mjl3hQj/2 I pj 



\mj\\ 



smh 

V 2 



1 -1 



\mj\ 



N- 



(24) 



for \mj\ > 1 and 



a^^ = -pJ coth i^^j ~ (25) 

for mj = 0. The density of modes 1]{Q) only contributes one additional factor of A^. So, 
if \mj\ > 1 the negative powers of the number of ions in the crystal cannot be compensated 
and the corresponding terms do not contribute. This result has already been used to derive 
equation ()14j) . For vanishing absolute temperature, the above expansion is always legitimate 
and one finds: 



a+i = +p'^ = -p-^ ai^ = 0. (26) 

Finally, the transition to the continuum has to be performed as described above. Now, 
equation (fTBj) can be reexpressed depending on the dimension d of the crystal: 



/,(0) = exp(-^ r\nn'-^ 

with the parameter: cr^ = {dqUY / {2MionhVL£,). The aforementioned separation into 
a 0-dependent and a (/(-independent part without leaving a non-integrable pole in the 0- 
dependent part is possible in three dimensions at any temperature and in two dimensions 



coth 



cos{0i7} — 1) — ism{(f)VL} 



(27) 
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for zero temperature. The analytic expression for the exponent Fa of the Mossbauer-Lamb 
fraction is given by: 



fdm-'-\oti.[^] (28) 



which for zero temperature is equal to F^{Q£)) = a^j id — V) and for non-zero temperature 
and in three dimensions reads: 



_ rin(l -e+^^^«) dilog(e+^^^^) + dilog(l-e+^^^«) vr^ 

where the dilogarithm is defined as in chapter 27.7 of reference^. For the subset of cases 
investigated up to the present, the exponent is bounded and the Fourier retransformation 
can be carried out term by term for the exponential series. The integrated contributions 
Id,u{4> = 0) follow a Poissonian distribution with a mean value equal to F^. So, this parameter 
determines how fast the series converges. For zero temperature the retransformations for all 
the terms of the series can be carried out exactly: 



(see figure H]) and 

(see figure 121) • The integrations in the three dimensional case and at finite temperatures 
can be carried out after approximating the hyperbolic cotangent by: 



coth 



^ (32) 



An expansion into partial fractions leads to (see figure : 



hA^) ^ infill Kx. (-^')'^-^e+^^^-'0(-c.') + 

Mi=0 M2=0 Ai=l 

+ E (W)^-^^(W) + 

A2 = l 

+ J2 ^3,A3 (W)^^"'e~'^^^«"'^(W) (33) 

A3 = l 
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with the factors independent of uj' = 2/i2 — /ii + u/Qd defined as: 



A'(l,2. 



_p 1 fpnnD 



(1,2,3), A 



and where the coefficients A(i^2,3),a are given by: 



A(1,2,3),A (-1)^ 
^^4.-Mi-A (A-1)! 



(34) 



Al,Ai 



1 



(Ai-1)! 



/,!/-Al 



(0) 



i(f)=+l3h 



A 



1 



^2(i/-^2)-A2 



2,A2 



A 



(A2 - 1)! d02(z.-^2)-A2 
1 rf^-'^^-Aa 



(_^0)2(^-M2)(_i0_^;i) 



U-fl2 



<j)=0 



-i(f) - (3hY-f'-^{+i(t) - (3hy 



3,A3 



(A3-I)! 

with the polynomial in 0: 



p 



(0) 



(35) 



p(0) = -z03(/5^f]^ + e-z^^f^^/^) - 02^^(1 - e-/3^f^^/2/2) - z0(/5^)2(/3^fi^/4) - (/3;i)V4 

(36) 

For the two-dimensional case at finite temperatures, the separation into a 0-dependent 
and a ^-independent part is no longer possible without leaving behind a non-integrable pole. 
Actually, the Mossbauer-Lamb fraction is exactly equal to zero because its exponent 
diverges. Hence, there is no recoil-free emission in this case. For further evaluation, we start 
by separating off the shift spectrum for zero temperature from that for finite temperature: 



0-2 



-I3hn 



cos{(f)Q} — r 



(37) 



With an approximation analogous to that in equation ()32|) one finds: 



,2 exp <j 2^^[E,{(3nnD/2) - Re{E,[{(3n/2 ± z0)r]^]}] 



where the exponential integral function Ei is defined as in chapter 5 of reference^: 



(38) 



Ei(z) 



dt- 



t 



(39) 



The following integrations cannot be performed analytically, due to the exponential in- 
tegral function Ei. For the two dimensional case we had already postulated that (T2 is to 
be small compared to one. In the physically interesting case of absolute temperatures T 
small compared to the Debye temperature 0_d = ^^D/ks, the exponent as such becomes 
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extremely small and the exponential function can be approximated by one. In the low- 
temperature approximation it remains to perform the Fourier retransformation of the power 
of the Lorentz distribution. However, the asymptotic form of this function for large is given 
by [/9/i/(20)]~^'^2/*^^^^^^. Due to the small absolute values of the exponent this function de- 
cays slower than which renders the Fourier retransformation difficult. Let us introduce 
an additional factor which ensures the convergence of the integral under investigation. If 
this manoeuvre is performed in a suitable way, the observable spectrum >S'(ci;^) can be ob- 
tained directly. When choosing a Lorentz distribution as a model for the natural line shape 
of the nucleus, the observabel spectrum is just the convolution of the former with the shift 
spectrum I{uj) in momentum representation. This corresponds to a multiplication of the 
Fourier transformed shift spectrum 7(0) with the Fourier transformed Lorentz distribution, 
i.e., a Laplace distribution also known as a doubly-exponential distribution. What remains 
to be executed is: 



°"2 



S.(u) ~ Pico) * / ^e^^(---)e-^l^l ^^"^ (40) 

with the resonance frequency of the nucleus ujq, the natural line-width ( (half-width half- 
maximum) and the integrated intensity 5*0. The integral can be interpreted as a Laplace 
transformation and yields (see-S' and figure Ej): 

^ * Re |l,_^ ( m^lt^) ] , (41) 



27r \ 2 0nsi 

with: 



U{z) = 2'-'V^T i^l + -J z~'[}li{z)-Yi{z)], (42) 

where B.i{z) denotes the Struve function (see chapter 12 of reference^), Yi{z) stands for a 
Bessel function of the second kind also known as Weber function (see chapter 9 of reference^), 
and r(^) represents the F-function. In the limit of vanishing line width (, the real part of 
the argument of l-ii{z) is equal to zero. Now, the approximation to what shall be called the 
Quasi- Mosshauer line of the shift spectrum is the convolution of the zeroth order term of 
the shift spectrum for zero temperature l2y=Q{oj) with the line shape given by equation (PT|) 
with 5*0, ojq, and C set to zero (see figure llj): 

^'""-^^ ' VLd \ \ 2 J } ^ ' 

IV. RESULTS 

In three dimensions there is largely agreement between the present work and-. Here, a 
deeper insight has been permitted due to the presentation of the results in an analytic form. 
The investigations have been extended to incorporate the two dimensional case at zero and 
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at finite temperature respectively. The one- dimensional case at zero temperature is also 
covered but not investigated any further here. 

In the two dimensional case at finite absolute temperatures, the Mossbauer-Lamb fraction 
according to its definition (|17|) vanishes exactly. No recoil-free absorption or emission takes 
place and strictly speaking, there is no Mossbauer line in the spectra. The mathematical 
investigation carried out here in this case is only applicable for small values of the model 
parameter ct2. But even then the execution of the calculations was not possible in the 
same way as for the other cases (three dimensions and/or zero absolute temperature). An 
alternative way had to be pursued. Although a Mossbauer line does not exist, a central peak 
is found. Already in the shift-spectrum representation this Quasi-Mossbauer peak shows a 
finite width. The width depends on the model parameter 02 and the temperature T. 

The disappearance of the Mossbauer-Lamb fraction in two dimensions and at finite ab- 
solute temperature signifies that no entirely recoil-free events are possible and that there 
cannot be a Mossbauer peak. Its replacement, the Quasi-Mossbauer line must be based 
on quasi-recoil-free processes. These are events during which only very little energy is 
transformed into lattice vibrations. This could be due to the predominant creation and an- 
nihilation of low-energy phonons. Another more general possibility is based on a dominance 
of multi-phonon events with small net energy. In the two-dimensional opposed to 

the three-dimensional, the density of modes does not drop rapidly enough for small absolute 
values of the circular frequency in order to compensate the pole of the shift amplitudes. 
This indicates that the first of the two explanations given above is the decisive one but 
does not rule out the second. The different mechanism of formation in two dimensions leads 
to a central peak of finite width depending on the temperature. In the representation as 
observable spectrum, i.e., after the convolution of the shift spectrum with the natural line 
shape, the Mossbauer and the Quasi-Mossbauer line appear qualitatively similar (see figure 
13}. Quantitatively they are not. They differ in peak height, line width, and integrated 
intensity. Given the qualitative similarity, the calculations would have to be reperformed 
with a more realistic than the Debye model in order to obtain more reliable quantitative 
results. 

While the phonon wings might render the determination of the central line's integrated 
intensity cumbersome, they provide the best means to identify two-dimensional elastic be- 
havior (see figure E)). For narrow spectra, in three dimensions, the phonon wing passes 
through a relative minimum almost reaching down to zero close to the central peak. In two 
dimensions a nearly horizontal passage can be seen. This means also that the phonon wings 
for positive and negative circular frequencies are continuously connected at zero frequency 
in three dimensions (see also figure El). In two dimensions there is a step (see figure Ej). 

Room for extension of the present work lies in the investigation of slowly decaying spectra 
in two dimensions. 



12 



Acknowledgements 

Help by R. J. Jelitto is gratefully acknowledged. 



G. K. Wertheim, Mossbauer Effect: Principles and Applications, Academic Press, New York 
(1964). 

2 H. J. Lipkin, Annals of Physics 9, 332-339 (1960). 

3 R. L. Mossbauer, Z. Naturforschg. 14 211-216 (1959). 
^ W. E. Lamb, Phys. Rev. 55 190 (1939). 

5 W. M. Visscher, Annals of Physics 9, 194-210 (1960). 

^ N. W. Ashcroft and N. D. Mermin, Solid State Physics, Chap. 23, Saunders College Publishing, 
(1976). 

^ M. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions, Dover, New York 
(1972). 

® A. Erdelyi, W. Magnus, F. Oberhettinger, and F. G. Tricomi, Tables of Integral Transforms, 

Chapter 4.3(10), McGraw-Hill, New York (1953). 
^ In thermal equilibrium: A{{n}i) = Y[j ■ 

^ The probability for the change of two configurations of the crystal into another is given by 
the product of the probabilities for the respective changes in all the modes: A^j^y{{n}^) = 
Uj a-', . 



13 



0.9 





0.6 - 




0.3 - 












1 



2 



3 



4 



5 



6 



7 



8 



9 



CO/Qd 



FIG. 1: The figure shows the shift spectrum I{uj) in three dimensions and for zero temperature. 
Different values of the model parameter have been chosen. The abscissa and the ordinate are 
rescaled with the Debye frequency Qd- In an unsealed presentation, a variation of the Debye 
frequency with the parameter kept constant would lead to variations of the aspect ratio of 
the graph. The step structure inherent to the shift spectrum with steps at integer multiples of the 
Debye frequency emerges from the plot. For smaller values of the Mossbauer-Lamb fraction these 
steps are smoothed out more and more until the spectrum resembles a Poissonian or Gaussian 
distribution. For large Mossbauer-Lamb fractions apart from the Mossbauer line which coincides 
with the axis in this plot, the spectrum mainly only consists of a small contribution between 
zero and the Debye frequency ^la- At zero temperature, the shift spectrum is identical to zero 
for negative values of the frequency u, because there is no initial population of the crystal with 
phonons which could be annihilated. Especially the sharp peak at lo = 0,^ can be traced back to 
the use of the Debye model and would be less pronounced for a model with an adapted behavior 
for higher phonon frequencies. 
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FIG. 2: Here, the shift spectrum I{io) for two dimensions is depicted for zero temperature, where 
various values have been chosen for the model parameter a2- The scaling properties and the step 
structure are equally well discernible as in figure ^ The transition to a Poissonian or Gaussian 
like distribution for slowly decaying spectra becomes evident, too. The difference to the three- 
dimensional case lies in the details of the different contributions. The spectrum for the highest 
Mossbauer-Lamb fraction in two dimensions shows an almost constant passage from zero frequency 
to the Debye frequency whereas in three dimensions it shows a linearly rising behaviour. 
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FIG. 3: In this figure, the shift spectrum I{io) is presented for different values of the absolute 
temperature T and the model parameter: = 2. For zero absolute temperature the corresponding 
plot is given in figure ^ The shift spectrum only depends on the ratio of the absolute temper- 
ature to the Debye temperature Qd- For the shown cases with relatively low temperatures, the 
step structure known from figure ^ is conserved. However, for increasing temperature it becomes 
smoother. The scaling behaviour with the Debye frequency VLd also remains the same, only that 
additionally, the ratio between the absolute and the Debye temperature Qd has to be kept con- 
stant. Now there are also contributions for negative circular frequencies w. They belong to those 
events where net energy is drawn from the crystal. 
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FIG. 4: The two figures 0] serve as a comparison between the Quasi-Mossbauer line in two dimen- 
sions and the Mossbauer line in three dimensions for the same choice of the physical constants: 
a"3 = 0.05 = (9/4)(T2, (3fi^D = 10, and C/^D = 10""^. The two graphs show the observable spec- 
trum Sobsi^^) rescaled with the Debye frequency J7d and the integrated intensity 5o. In order to 
emphasise the differences, a small natural line-width has been chosen. The lines are always even 
functions with respect to inversion of the sign of the argument relative to the resonance frequency 
loq. The figure on the Ihs is a doubly linear plot. Here, the slightly higher maximum value of 
the Mossbauer line (dashed) compared to the Quasi-Mossbauer line (solid) can be seen. For large 
values of the circular frequency, the Mossbauer line tends to zero faster than the Quasi-Mossbauer 
line. This can best be seen in the figure on the rhs which has a logarithmic ordinate. The two 
lines do not have the same integrated area. The area beneath the Mossbauer line is larger by a 
factor of e^'^/^^'^^. On first sight, the two lines do not differ in their width. Numerically a difference 
of a Quasi-Mossbauer line with a 0.4 percent larger width is found. This value can be increased 
to 3 percent for the smallest natural line width found in nature {C/Qd = 10~^^). In the last case 
however, the absolute difference is one million times smaller. 
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FIG. 5: In this figure, the graphs of the following two observable spectra Sobsi^) at finite tem- 
perature are depicted: The one for two dimensions (dashed) with the model parameter set to 
(72 = 1/45 and the other for three dimensions (solid) with the model parameter = 1/20. The 
other parameters defining the temperature and the natural line-width respectively are: phfto = 10 
and C/^D = 10~^. The spectrum in two dimensions decays more slowly for negative values of the 
circular frequency u than in three dimensions. The most pronounced difference can be found for 
the circular frequency between zero and the Debye frequency Qd- There, the contribution of the 
first order of the expansion of the exponential function shows up. In three dimensions, the ob- 
servable spectrum Sobsi^) passes through a relative minimum almost reaching down to zero. In 
two dimensions, a nearly horizontal passage can be seen. This interval is the most suitable place 
to distinguish between the two cases by measurements and to identify two-dimensional elastic 
behaviour. 
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